function Qn = fn_u2Qn(n, H, T, un)

a = fn_inveta(n, H/T);
Qn = a./(exp(un) + a);

end